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We consider a quantum impurity model in which a bosonic impurity level is coupled to a non- 
interacting bosonic bath, with the bosons at the impurity site subject to a local Coulomb repulsion 
U . Numerical renormalization group calculations for this bosonic single-impurity Anderson model 
reveal a zero-temperature phase diagram where Mott phases with reduced charge fluctuations are 
separated from a Bose-Einstein condensed phase by lines of quantum critical points. We discuss 
possible realizations of this model, such as atomic quantum dots in optical lattices. Furthermore, 
the bosonic single-impurity Anderson model appears as an effective impurity model in a dynamical 
mean-field theory of the Bose-Hubbard model. 

PACS numbers: 05.10.Cc (Renormalization Group methods), 05.30.Jp (Boson systems), 03.75.Nt (Other 
Bose-Einstein condensation phenomena), 03.75.Hh (Static properties of condensates). 



The focus of this work is the physics of a bosonic im- 
purity state coupled to a non-interacting bosonic envi- 
ronment modeled by the Hamiltonian 

H = eQb^b+^Ub^b{b^b-l) 

+ J2ekbibk + J2Vk{blb + b^bk) . (1) 

k k 

The energy of the impurity level (with operators 6^1^) is 
given by £o; the parameter U is the local Coulomb re- 
pulsion acting on the bosons at the impurity site. The 
impurity couples to a bosonic bath via the hybridization 
Vk, with the bath degrees of freedom given by the oper- 
ators b'^iJ'^ with energy e^. 

We term the system defined by eq. the 'bosonic 
single-impurity Anderson model' (bosonic siAm) , in anal- 
ogy to the standard (fermionic) siAm [3] which has a very 
similar structure except that all fermionic operators are 
replaced by bosonic ones. Furthermore, we do not con- 
sider internal degrees of freedom of the bosons, such as 
the spin (an essential ingredient in the fermionic siAm). 

There are various ways to motivate the study of the 
model eq. From a purely theoretical point of view 
it is interesting to compare the physics of the fermionic 
and bosonic versions of the siAm. Of course, there are 
striking differences between these two models: there is 
no direct bosonic analog of local moment formation and 
screening of these local moments at low temperatures, 
characteristics of the fermionic Kondo effect 0- On 
the other hand, the bosonic model allows for a Bose- 
Einstein condensation (BEG), at least in a certain pa- 
rameter space (see Fig. Q] below), a phenomenon which is 
clearly absent in the fermionic model. 

There are certain similarities between the fermionic 
and the bosonic model concerning the role of the local 
Coulomb repulsion U: increasing the value of U can 
induce a quantum phase transition from a phase with 
screened spin to a local moment phase in the fermionic 
case while it can induce a quantum phase transition 



from a BEC phase to a 'Mott phase' in the bosonic case, 
as shown below. 

Another motivation for studying the bosonic siAm 
comes from a treatment of the Bose-Hubbard model 
within dynamical mean-field theory (DMFT) 01 ■ Al- 
though such an investigation has not been pursued so 
far, it is clear that the effective impurity model onto 
which the Bose-Hubbard model is mapped will have a 
similar form as eq. (see also Ref. 0) . An obvious ap- 
plication of such a DMFT treatment would then address 
the Mott transition in the Bose-Hubbard model which 
must have its counterpart on the level of the effective im- 
purity model - and the 'Mott transition' on the impurity 
level is precisely what we are looking at here. 

Finally, a physical system described by the bosonic 
siAm could be directly realized in optical lattices, where 
the laser fields are tuned in such a way that a single 
'impurity' site is formed within an approximately un- 
perturbed lattice system ('atomic quantum dots', see 
Refs.0,0)- Provided the Coulomb repulsion at all sites 
except the impurity can be neglected, the corresponding 
model can be directly mapped onto the model eq. Q). 
Present theoretical studies of atomic quantum dots focus, 
however, on a coupling between impurity and excitations 
of the superfluid environment , for which a description 
in terms of the spin-boson model is more appropriate. 

For the calculations presented in this paper we use the 
numerical renormalization group (NRG) originally devel- 
oped by Wilson for the Kondo problem Q . This method 
has been shown to give very accurate results for a broad 
range of impurity models, including the fermionic siAm 
P, 0| but also quantum impurities with coupling to a 
bosonic bath 0, ^| . Here we employ the bosonic ex- 
tension of the NRG (bosonic NRG) to study the model 
eq. (QJ. We shall present ground state phase diagrams 
of the bosonic siAm, identify fixed points and discuss 
quantum phase transitions between the fixed points. Al- 
though far from comprehensive, our results indicate that 
the bosonic siAm shows a variety of interesting proper- 



ties which deserve to be studied in greater detail in the 
future. 

Before we come to the results from the bosonic NRG, 
let us discuss some general properties and trivial limits 
of the model eq. (Q. 

Similar to other quantum impurity models, the influ- 
ence of the bath on the impurity is completely specified 
by the bath spectral function 



(2) 



Here we assume that A(w) can be parametrized by a 
powerlaw for frequencies up to a cutoff lUc (we set lUc — I 
in the calculations) 



A{uj) ^ 2t: a LO^ ^ , < w < w, 



(3) 



The parameter a is the dimensionless coupling constant 
for the impurity-bath interaction. This form of A.(uj) is 
certainly not the most general one and specific applica- 
tions (within a bosonic DMFT or for an impurity level 
in an optical lattice) will lead to additional structures in 
A(c^). 

The NRG calculations can be both performed for 
a canonical and a grandcanonical ensemble. Here we 
present only results for a grandcanonical ensemble for 
which we set the chemical potential of the bath to /i = 0. 
This means that a Bose-Einstein condensation of the 
whole system (impurity plus bath) can only be induced 
by the coupling to the impurity. This can already be 
understood from the non-interacting case, U = 0: here 
a direct diagonalization of the bosonic siAm shows that 
with increasing a, a localized state with negative energy 
separates out of the continuum at a critical value a = ac- 
The BEC then occurs via populating this localized state, 
a feature which can also be observed in the numerical 
calculations. 

The other trivial limit of the bosonic siAm is the decou- 
pled impurity, a = 0. The succession of quantum phase 
transitions as shown in Fig. ^ for a — can be easily un- 
derstood from the dependence of the many-particle levels 
on the parameters Eq and U. The transition occurs for 
£a/U — — riinip when the energies of the states with rtimp 
and nimp + 1 electrons are degenerate. 

The full phase diagram Fig. ^ is calculated with the 
bosonic NRG In this approach, the frequency range 
of the bath spectral function [0,a;c] is divided into in- 
tervals [wcA-("+i\wcA~"], n = 0,1,2,..., with A the 
NRG discretization parameter (we use A = 2.0 for all 
the results shown in this paper). The continuous spec- 
tral function within these intervals is approximated by a 
single bosonic state and the resulting discretized model is 
then mapped onto a semi-infinite chain with the Hamil- 
tonian 




FIG. 1: Zero-temperature phase diagram of the bosonic siAm 
for bath exponent s = 0.6 and fixed impurity Coulomb in- 
teraction U = 0.5. The different symbols denote the phase 
boundaries between the Mott phases and the BEC phase. The 
Mott phases are labeled by their occupation nimp for a = 0. 
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parameters are A — 2.0, A'^b ~ 10, 
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Here the impurity couples to the first site of the chain via 
the hybridization V = ^/2a/(T+s). The bath degrees 
of freedom are in the form of a tight-binding chain with 
operators bn\ on-site energies e„, and hopping matrix 
elements i„ which both fall off exponentially: tn,en oc 
A~". 

The chain model eq. Q is diagonalized iteratively 
starting with the impurity site and adding one site of 
the chain in each iteration. As for the application to the 
spin-boson model, only a finite number of basis states 
for the added site can be taken into account, and, after 
diagonalizing the enlarged cluster, only the lowest-lying 
Ng many-particle states are kept for the subsequent iter- 
ations (here we use A^^b = 10 - 20 and = 100 - 200). 
The main technical difference to the spin-boson model 
is that we can use the total particle-number as a con- 
served quantity in the Hamiltonian eq. Furthermore, 
the renormalization group flow turns out to be consider- 
ably more stable as in the calculations for the spin-boson 
model so that we can easily perform up to A^ = 100 it- 
erations (a detailed account of the technical details will 
appear elsewhere). 

Let us now discuss the T — phase diagram of the 
bosonic siAm, Fig. ^ calculated for fixed U — 0.5 with 
the parameter space spanned by the dimensionless cou- 
pling constant a and the impurity energy Sq. We choose 
s — 0.6 as the exponent of the powerlaw in A(tt') (the s- 
dependence of the phase diagram is discussed in Fig.0]be- 
low). The phase diagram is characterized by a sequence 
of lobes which we label by the impurity occupation rtimp 
at a = where it takes integer values. We use the termi- 
nology 'Mott phases' for these lobes, due to the appar- 
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ent similarity to the phase diagram of the Bose-Hubbard 
model. The Mott phases are separated from the BEC 
phase by lines of quantum critical points which termi- 
nate (for s — 0.6) at a finite value of a, except for the 
n-imp = phase where the boundary extends up to infi- 
nite a. These transition can be viewed as the impurity 
analogue of the Mott transition in the lattice model, since 
it is the local Coulomb repulsion which prevents the for- 
mation of the BEC state. 

The phase diagram Fig. ^ is deduced from the flow of 
the lowest-lying many-particle levels which allows quite 
generally to identi fy fi xed points and transitions between 
these fixed points [lOj. Figure |5] shows the fiow for s = 
0.4, fixed a = 0.007 and U = 0.5, and two values of 
So very close to the quantum phase transition. The solid 
lines belong to the Mott phase with ni,„p = 2 and display 
a crossover between two fixed points: from an unstable 
quantum critical point for iteration numbers N ~ 20 — 40 
to a stable fixed point for > 70. Analysis of the stable 
fixed point shows that it can be described by a decoupled 
impurity with occupation riimp = 2 and a free bosonic 
bath given by the decoupled chain. (The actual impurity 
occupation, however, differs from the integer values, as 
discussed below.) The structure of the quantum critical 
point, on the other hand, is presently not clear but might 
be accessible to perturbative methods as discussed in [l^ . 

The dashed lines belong to the BEC phase where, ap- 
parently, something dramatic is happening at iV « 55. 
The divergence of the energies for all excited states seen 
in this plot is due to the formation of a localized state, 
split off from the continuum by an energy gap Ag, very 
similar to the behavior for U ~ 0. This energy gap takes 
a finite value which is not renormalized to zero as 
as the levels of the other fixed points, therefore the di- 
vergence of Ag • A^. In the inset of Fig. |2] we plot En 
(instead of E^A^) so that the development of the gap 
Ag ?a 4 X 10~^^ appears as a plateau of the first excitated 
state. The extremely small value of the gap is due to the 
tuning of the parameters very close to the transition (Ag 
vanishes at the transition). 

Further analysis of the data shows that the ground 
state in the BEC phase has an occupation number given 
by the maximum boson number used in the iterative di- 
agonalization. 

The fact that the unstable fixed point separating the 
Mott phase from the BEC phase is indeed a quantum 
critical point can be deduced from its non-trivial level 
structure, as mentioned above, and the behavior of the 
crossover scale T*. Numerically we find that upon vari- 
ation of So close to its critical value eo,c, the crossover 
scale vanishes with a powerlaw at the transition, T* oc 
ko ~ £o,c|'^ on both sides of the transition, with a non- 
trivial exponent v w 2.50 ± 0.06. This exponent is ob- 
served for all lines of quantum critical points for fixed 
s = 0.4 and preliminary results show that — l/s for 
< s < 1. 




FIG. 2: Flow diagram of the lowest lying many-particle levels 
Em versus iteration number N for parameters s = 0.4, a = 
0.007, U — 0.5, and two values of eo very close to the quantum 
phase transition between the Mott phase with nimp ~ 2 and 
the BEC phase. Both the quantum critical point and the 
Mott phase appear as fixed points in this scheme whereas in 
the BEC phase, a gap Ag opens between the ground state 
and the first excited state, see the inset where En (instead of 
EpfA.'^) is plotted versus A''. 

Let us now focus on the impurity occupation Tiimp for 
temperature T = 0. Figure 01 shows the dependence of 
nimp on £o for s = 0.4, U = 0.5, and various values of a. 
The symbols indicate that the parameters lie within the 
Mott phase whereas the crosses are for the BEC phase. 
Both sets of lines form a continuous curve but it is unclear 
at the moment whether the data in the BEC phases are 
reliable (due to the special features in the flow diagram 
discussed above). 

The terminology we use here suggests an integer oc- 
cupation throughout the Mott phase, as for the Bose- 
Hubbard model, but for the bosonic siAm nimp deviates 
from the integer values as soon as the coupling to the 
bath is finite, see Fig. |3| This is to be expected since for 
a gapless bath spectral function A(ci;), the charge fluc- 
tuations on the impurity site cannot be completely sup- 
pressed. Indeed we observe that for increasing the value 
of the bath exponent s, the nimp(£o)-curve gets closer to 
the step function. At this point one can speculate about 
the possible development of A{uj) in a DMFT treatment 
of the Mott phase. The self-consistency might generate 
a bath spectral function with a gap and the impurity oc- 
cupation might then turn into the step function expected 
for the lattice model. 

The precise shape of the boundaries in the phase di- 
agram Fig. n] depends on the form of A{uj) for all fre- 
quencies. Here we stick to the powerlaw form eq. (O and 
present the dependence of the phase diagram on the bath 
exponent s in Fig.0] We observe that upon increasing the 
value of s, the areas occupied by the Mott phases extend 
to larger values of a and significantly change their shape. 
A qualitative change is observed when the exponent ap- 
proaches s = 1. First of all, the Mott phases appear 



4 



2 



1.5 



0.5 





a=0 

G — G 0=0.007 
■ ■ 0=0.010 




J 




4i 


» « 0=0.028 




/ . 














+" » j 

* :\ 

e — »i ■ ^ ' ' 









-05 0.5 1 1.5 2 

-£„/U 

FIG. 3: Impurity occupation function of eo for 

temperature T = 0, s = 0.4, U — 0.5, and various values of a. 
The sharp steps for the decoupled impurity a = are rounded 
for any finite a. Symbols (crosses) correspond to data points 
within the Mott (BEG) phases. 
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FIG. 4: Zero-temperature phase diagram of the bosonic si Am 
as in Fig. Q but now for different values of the bath expo- 
nent s. For increasing value of s, the areas occupied by the 
Mott phases significantly change their shape and for s = 1 it 
appears that each Mott phase extends up to arbitrarily large 
values of a. 

to extend up to arbitrarily large values of a. Further- 
more, the BEC phase which separates the Mott phases 
for s < 1 and a > 1 is completely absent for s = l! We 
do not yet have an explanation for this observation and 
it would be interesting to find out whether the absence of 
the BEC phase is due to the special form of A(a;), eq. l|3Jl, 
or whether it is a generic feature even when Aiuj) oc w is 
only valid for w — > 0. 

The case s = (constant bath density of states) turns 
out to be difficult to access numerically. An extrapola- 
tion of the phase boundaries for values of s in the range 
0.1...0.4tos = 0is inconclusive, but the Mott phase is 
at least significantly suppressed in this limit. 

To summarize, we have presented NRG calculations 
for the phase diagram and the impurity occupation of a 
bosonic version of the single-impurity Anderson model. 
The phase diagram contains Mott phases, in which the 
local Coulomb repulsion prevents Bose-Einstein conden- 



sation, separated by the BEC phase by lines of quan- 
tum critical points. The studies presented here are only 
the starting point for a comprehensive investigation of 
the bosonic siAm and we are planning to calculate, for 
example, physical properties at finite temperatures and 
dynamic quantities (impurity spectral function and self- 
energy). The latter will be of importance for a possible 
DMFT for the Bose-Hubbard model, an approach which 
has not yet been fully developed due to various concep- 
tual problems. One issue is the proper scaling of the 
Hamiltonian parameters in the limit of infinite spatial 
dimensions For example, the model on a hypercu- 
bic lattice as studied in Ref. |la requires a scaling of the 
hopping matrix elements as 1/d which leads to a static 
mean-field theory. In addition, the bosonic DMFT in the 
superfluid phase of the Bose-Hubbard model might gen- 
erate a more complex impurity model, the bosonic siAm 
introduced here would then be applicable only within the 
Mott phases of the lattice model. 

Finally, it would be interesting to identify situations for 
atomic quantum dots in optical lattices which can be de- 
scribed by the bosonic single-impurity Anderson model. 
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